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Abstract 

This paper deals with conservation laws on networks, represented 
by graphs. Entropy-type conditions are considered to determine dy- 
namics at nodes. Since entropy dispersion is a local concept, we con- 
sider a network composed by a single node J with n incoming and 
m outgoing arcs. We extend at J the classical Kruzkov entropy ob- 
taining two conditions, denoted by (El) and (E2): the first requiring 
entropy condition for all Kruzkov entropies, the second only for the 
value corresponding to sonic point. First we show that in case n 7^ m, 
no Riemann solver can satisfy the strongest condition. Then we char- 
acterize all the Riemann solvers at J satisfying the strongest condition 
(El), in the case of nodes with at most two incoming and two outgo- 
ing arcs. Finally we focus three different Riemann solvers, introduced 
in previous papers. In particular, we show that the Riemann solver 
introduced for data networks is the only one always satisfying (E2). 

Key Words: scalar conservation laws, traffic flow, Riemann solver, net- 
works, entropy conditions. 

AMS Subject Classifications: 90B20, 35L65. 

*E-mail: mauro.garavello@mfn.unipmn.it. Partially supported by Dipartimento di 
Matematica e Applicazioni, Universita di Milano-Bicocca. 
^E-mail: piccoli@iac.rm.cnr.it. 



1 Introduction 



Nonlinear hyperbolic conservation laws on networks have recently attracted 
a lot of interest in various fields: car traffic P, [161 lEl 122] , gas dynamics [U 
[21[31[71[81E1[I01[II1^[I3], irrigation channels [2 [El [201 El] and supply 
chains [SI [21] ■ A network is modeled by a graph: a finite collection of arcs 
connected together by vertices. On each arc we consider a scalar conservation 
law. For instance one may think to the Lighthill-Whitham-Richards model 
for car traffic [25l[2Z]- However, our results applies to the other application 
domains. 

It is easy to check that the dynamic at nodes is not uniquely determined 
by imposing the conservation of mass through vertices. Then, to completely 
describe the network load evolution, the first step is to appropriately define 
the concept of solution at a vertex. 

As in the classical theory of conservation laws, this problem is equivalent 
to giving the solution Riemann problems (now at vertices). More precisely, 
a Riemann problem at a vertex is simply a Cauchy problem with constant 
initial conditions in each arc of the vertex. The map, which associates the 
solution to each Riemann problem at a vertex J, is called a Riemann solver 
at J. Similarly to the case of a real line, one has to resort to the concept 
of weak solutions in the sense of distributions and there are infinitely many 
Riemann solvers producing weak solutions. First one uses entropy type con- 
ditions inside arcs as for the real line. Then, in order to select a particular 
solution (i.e. a Riemann solver) at the vertex, one has to impose some addi- 
tional conditions. In [6], for example, the authors required some rules about 
the distribution of the fiuxes in the arcs and a maximization condition; see 
also [2S] ■ It is then natural to ask if entropy-like conditions can be im- 
posed also at the vertex and not only inside arcs. 

In this paper, we focus on a single vertex J, composed by n incoming and 
m outgoing arcs and we extend the Kruzkov [23j entropy-type conditions. 
More precisely, we propose two different entropy conditions for admissibil- 
ity of solutions, called, respectively, (El) and (E2). The condition (El) is 
stronger than (E2), indeed the first asks for Kruzkov entropy condition to be 
verified for all entropies, while the second asks only for the precise Kruzkov 
entropy corresponding to sonic point. It is interesting to note that the en- 
tropy condition (El) imposes strong restrictions both on Riemann solvers 
and on the geometry of the vertex. Indeed, Riemann solvers satisfying (El) 
can exist only in the case of vertices with the same number of incoming and 
outgoing arcs. 

We then test our conditions on Riemann solvers considered in the litera- 
ture. First we can prove that the Riemann solver, introduced in [H] for data 
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networks, satisfies (E2) and, in special situations, also (El). 
Then we show that the Riemann solvers defined in [HIES] do not satisfy (E2). 
However, at least for the Riemann solver in [6], the entropy condition and 
the maximization procedure agree on some particular set, over which the 
maximization is taken. Roughly speaking the solver respects the entropy 
condition once traffic distribution is imposed. 

The paper is organized as follows. Section [2] introduces the basic defini- 
tions of networks and of solutions. Section [3] deals with the solution to the 
Riemann problem at the vertex J. Moreover, we introduce the entropy con- 
ditions (El) and (E2) for Riemann solvers at J. In Section HI we determine 
which Riemann solvers satisfy the entropy condition (El). The paper ends 
with Sectional which considers the Riemann solvers TZSi, IZS2 and IZS^, 
introduced respectively in [51 [HI [2S] , and analyzes what entropy conditions 
these Riemann solvers satisfy. 

2 Basic Definitions and Notations 

Consider a node J with n incoming arcs and m outgoing arcs 

J„+i, . . . , In+m- We model each incoming arc Jj (i G {1, . . . , n}) of the node 
with the real interval Jj =] — 00, 0] and each outgoing arc Ij (j G {n + 
1, . . . , n + m}) of the node with the real interval Ij = [0, +oo[. On each arc 
/; (/ G {1, . . . , n + m}), the traffic evolution is given by 

{pi)t + fiPi). = 0, (1) 

where pi = pi(t, x) G [0, Pmax], is the density, vi = Vi{pi) is the average velocity 
and f{pi) = vi{pi) pi is the flux. Hence the network load is described by a 
finite collection of functions pi defined on [0, +oo[xli. For simplicity, we put 
Pmax = 1- On the fiux / we make the following assumption 

(JF) / : [0, 1] ^ M is a piecewise smooth concave function satisfying 

1. /(0) = /(l) = 0; 

2. there exists a unique cr g]0, 1[ such that / is strictly increasing in 
[0, a[ and strictly decreasing in ]a, 1]. 

Definition 2.1 Let r : [0, 1] [0, 1] be the map such that: 
1- fi^ip)) = fip) for every p G [0, 1]; 
2. t{p) 7^ p for every p G [0, 1] \ {a}. 
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Definition 2.2 A function pi G C([0, +oo[; L]^^{Ii)) is an entropy- admissible 
solution to (QP in the arc Ii if, for every k G [0, Pmax] and every (p : [0, +oo[xli 
R smooth, positive with compact support in ]0, +oo[x (/; \ {0}) 

j^^ j^^ {\pi~k\^ + sgn(p, - k)U{pi) - fm^)dxdt > 0. (2) 

Definition 2.3 A collection of functions pi G C{\Q,+oo[;L]^J^Ii)), (I G 
{1, . . . , n + m} ) is a weak solution at J if 

1. for every I G {1, . . . ,n + m}, the function pi is an entropy- admissible 
solution to (QP in the arc Ii; 

2. for every / G {1, . . . , n-\-m} and for a.e. t > 0, the function x i-^ pi{t, x) 
has a version with bounded total variation; 

3. for a.e. t > 0, it holds 

n n+m 

E/(/^^(^'0-))= E M(^'0+))' (3) 

i=l j=n+l 

where pi stands for the version with bounded total variation. 



3 The Riemann Problem at J 

Given pi_o, • • • , Pn+mfi ^ [0, 1], a Riemann problem at J is a Cauchy problem 
at J with constant initial data on each arc, i.e. 

I e {l,...,n + m}. (4) 

P«(0, ■) = Po,h 

Now, we give some definitions for later use. The first one is the defini- 
tion of Riemann solver, which is a map giving a solution to the Riemann 
problem (jl]). 

Definition 3.1 A Riemann solver TZS is a function 

lis : [0, 1]'^+"^ — > [0, 1]''+'" 

(Pl,0! • • • ) Pn+m,o) ' ^ (Pl) • • • ) Pn+m) 

satisfying 



2. for every i G {1, . . . , n}, the classical Riemann problem 
{ pt + fip)x = 0, xGM, t>0, 



p(0,x) = | 




is solved with waves with negative speed; 
3. for every j G {n + 1, . . . ,n + m}, the classical Riemann problem 



is solved with waves with positive speed. 

We introduce the concepts of equilibrium and consistency for Riemann 
solvers. The fixed points of a Riemann solver are called equilibria, while a 
Riemann solver has the consistency condition when its image is contained in 
the equilibria. 

Definition 3.2 We say that (pi,o, ■ • • , Pn+m,o) is an equilibrium for the Rie- 
mann solver TZS if 



Definition 3.3 We say that a Riemann solver IZS satisfies the consistency 
condition if, for every (pi,o, . . . , p„+m,o) e [0, 1]"+'", then 7^5(pl,o, . . . , Pn+m,o) 
is an equilibrium for TZS . 

We introduce now the concepts of entropy functions and admissible en- 
tropy conditions (El) and (E2) for Riemann solvers. We are essentially 
extending the Kruzkov entropy condition to the case of a node; see |23] . 



Definition 3.4 The function T : [0, 1]'^+'^ x [0, 1] M, defined by 




x G M, t > 0, 



) = (Pl 



n 





i=l 



n+m 



J2 sgn(p, - k) (/(p,) - f{k)) , 



j=n+l 



is called entropy-flux function. 
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Definition 3.5 A Riemann solver TZS satisfies the entropy condition (El) 
if, for every initial condition (pi,05 • • • ,Pn+m,o) Oj'i^d for every k G [0, 1], we 
have 

T{pi, . . .,Pn+m, k) > 0, (6) 
where (pi, . . . , pn+m) = 7?.5(pi,o, • • • , p„+m,o)- 

Remark 1 If k = 0, then equation ^ becomes Yl'^=i fiPi) — Yl'j=n+i fiPj)- 
If k = 1, then equation ^ becomes Yll=ifiPi) — ^]=n+i f(Pj)- Therefore 
the entropy condition (El) implies the conservation identity J2i=ifi.Pi) = 
EJ=r+i/(Pi)- 

Definition 3.6 A Riemann solver TZS satisfies the entropy condition (E2) 
if, for every initial condition (pi,o, • • • ,Pn+m,o); we have 

-^(Pi, • • • ,Pn+m,o-) > 0, (7) 

where (pi, . . . , pn+m) = TZS{pi^o, . . . , p„+m,o)- 

Remark 2 The entropy condition ^ can be deduced in the following way. 

Fix, for every I G {1, . . . , n + m}, a smooth function (fi : [0, +oo[x// — > 
[0, +oo[ with support contained in [0, +oo[x [— M, M] for some M > and 
assume that ipi'{t, 0) = fi"(t, 0) for every t > and I', I" G {1, . . . , n + m}. 
Applying the divergence theorem to the inequality 

n+m „+oo p 

y2 / [\Pi-k\'~pi,t + sgn{pi-k){f{pi)- f{k))ifijdxdt>0, 

i^i Jo Jli 

where (pi, . . . ,pn+m) is an equilibrium at J, we deduce 

Obviously, these kinds of entropies are not justified by physical consider- 
ations. 

Finally, let us introduce sets fli and related to the points [2] and [3] of 
Definition I3.1[ 

1. For every i G {1, . . . ,n} define 

[0, /(Pi,o)], if < Pi,o < cr, 
[0,/((t)], ifa<p,,o<l, ^ 



and 



, {pi,o}U]r(pi,o), 1], if < pi,o < CT, 
[0^,1], if 0- < Pi,0 < 1- 
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2. For every j G {n + 1, . . . ,n + m} define 



^^_^|o./W]. ifo<p,,„<.. ^^^^ 

[0,/(Pj,o)], if 0^ < Pj,o < 1, 



and 



, [0,c^], ifO<Pj,o<fT, 
<l>j = <^ (11) 
{pjfl} U [0, r(pj-o)[, if (T < pj-o < 1. 

Tlie following Proposition links the previous sets with Definition 13.11 

Proposition 3.1 The following statements hold. 

1. For every i G {1, . . . ,n} , an element 7 belongs to Qi if and only if there 
exists pi G [0, 1] such that f{pi) = 7 and pointlE of Definition \3.1\ is 
satisfied. 

2. For every j G {ra + 1, . . . , n + m}, an element 7 belongs to Qj if and 
only if there exists pj G [0, 1] such that f{pj) = 7 and point of 
Definition \3.1\ is satisfied. 

The proof is trivial and hence omitted. The main result of this Section 
is that, if n 7^ m, then every Riemann solver TZS at J does not satisfy the 
entropy condition (El). We first need the following result. 

Proposition 3.2 Fix a node J with n incoming arcs and m outgoing arcs 
and a Riemann solver TZS satisfying the entropy condition (El). Denote with 
(pi, . . . ,p„+m) the image throughlZS of the initial condition (pi,o, • • • ,Pn+m,o)- 

1. If n > m, then mm {pi, pn} = 0. 

2. Ifn< m, then max {p„+i, . . . , Pn+m} = 1- 

Proof. Consider first the case n > m. Suppose by contradiction that 
min {pi, . . . , pn} > 0. Define the set J = {j E {n + 1, . . . , n + m} : pj = 0} 
and fix < A; < min {pi : / G {1, . . . , n + m} \ J}. Thus, the entropy in- 
equality JF (pi, . . . , pn+m, k) > becomes, 

n 

E [/(/^^) - /(^)] ^ E if(p^^ - /(^)] + E /(^)- 

i=l j€{n+l,...,n+m}\.J j£J 

By point [H of Definition 13.11 we deduce that 

-nf{k) > -{m - MJ))f{k) + #(J)/(fc), 
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where #(^7) denotes the cardinahty of J; thus {m — n — 2#(J))/(A;) > 0, 
which is a contradiction. 

Consider now the situation n < m. By contradiction we assume that 
max {pn+i, • • • , pn+m} < 1- Define the set / = {i G {1, . . . , n} : pi = 1} and 
fix max {pi : I E {I, . . . ,n + m} \ 1} < k < 1. Thus, the entropy inequahty 
JT (pi, . . . , pn+m, k) >0 becomes, 

n+m 

je{l,...,n}\/ ie/ j=n+l 

By point [T]of Definition 13. 11 we deduce that {n — 2^{I) — m)f{k) > 0, which 
is a contradiction. □ 



Theorem 3.1 Fix a node J with n incoming arcs and m outgoing arcs and 
suppose that n ^ m. Every Riemann solver TZS at J does not satisfy the 
entropy condition (El). 

Proof. Suppose, by contradiction, that there exists a Riemann solver TZS 
at J satisfying the entropy condition (El). 

Assume n> m and consider an initial condition (pi,o, • • • , Pn+m,o) satisfy- 
ing pi^o ^ for every z G {1, . . . , ra}. If (pi, . . . , pn+m) = 7^5(pi,o, • • • , Pn+m,o), 
then, by Proposition 13.21 there exists zi G {1, . . . , n} such that pi^ = 0, which 
is a contradiction since the wave {pii,o,Pii) has not negative speed. 

Assume now n < m and consider an initial condition (pi,o, • • • ,Pn+m,o) 
satisfying pj Q ^ 1 for every j G {n + 1, . . . , n + m}. By Proposition 13.21 if 
(pi, . . . , pn+m) = 7^5(pi,o, . . . , Pn+m,o), then there exists ji e {n + 1, . . . ,n + 
m} such that pj^ = 1, which is a contradiction since the wave (pji,Pji,o) has 
not positive speed. □ 



4 Riemann solvers satisfying (El) 

In this Section we determine which Riemann solver satisfies the entropy con- 
dition (El), in the sense of Definition 13.51 for nodes with n = m G {1,2}. 
In the case n ^ m, Theorem 13.11 implies that every Riemann solver does not 
satisfy (El). Moreover if n = m = 1, then there exists exactly one Riemann 
solver at J satisfying (El), while if n = m = 2, then there exist infinitely 
many Riemann solvers satisfying (El); see Sections 14.11 and 14.21 We do not 
treat the case n = m > 2, for the huge number of different situations. 
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4.1 Nodes with n = m = 1 



In this subsection, we fix a node J with one incoming and one outgoing arc. 
The following result holds. 

Proposition 4.1 A Riemann solver TZS at J satisfies the entropy condition 
(El) if and only if, for every initial datum (pi,o, P2,o); th^ image (pi,p2) = 
7lS{pifi, p2fl) satisfies either 

Pi=P2 (12) 

or 

Pi < p2 and f{pi) = f{p2). (13) 

Proof. Consider first a Riemann solver TZS satisfying the entropy condition 
(El). By [U of Definition 13.11 it is clear that /(pi) = /(p2)- Assume by 
contradiction that pi > p2- Since /(pi) = /(P2), we easily deduce that 
P2 < o" < pi. Putting A; = cr in equation ([6]) we derive 

/(Pl)-/(a)>/(a)-/(p2), 

which is, by assumptions, equivalent to /(pi) > /(cr), and so we get a con- 
tradiction. 

Consider now a Riemann solver TZS such that, for every initial datum 
(pi,o,P2,o), the image (pi,P2) = 7?.5(pi,o, P2,o) satisfies either or ([I3l). It 
is trivial to prove that (El) holds. □ 



Theorem 4.1 There exists a unique Riemann solver TZS at J satisfying 
the entropy condition (El). This Riemann solver satisfies the consistency 
condition and coincides with the Riemann solver introduced in for traffic 
or with the Riemann solver introduced in |7^. 

Proof. Fix an initial datum (pi,05P2,o)- We show that there exists a unique 
(pi,P2), which is the image of an entropy admissible Riemann solver. 

If pifi = p2,o, then we claim that pi = P2 = Pi,o- Assume by contradiction 
that pi 7^ p2- In this case either pi < cr < p2 or p2 < cr < pi. By Propo- 
sition 14.11 the only possibility is pi < cr < p2. By Proposition 13.11 either 
Pi = Pi,o or P2 = P2,o- In the first case p2 = t(p2,o); while in the second one 
Pi = t(pi,o)- It is not possible. 

Assume now that pi^ 7^ p2,o- We have some different possibilities. 

1. max{pi 0, P2,o} ^ o"- By Proposition 13. H we deduce that p2 G [0,a]. 
Moreover, by Proposition 14.11 we deduce that pi = pi^o! hence p2 = 
Pi = Pi,o- This solution respects all the properties of Definition 13. II and 
the entropy condition ([6]). 
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2. min{pi 0, P2,o} > cr- By Proposition \3.1\ we deduce that pi G [cr, 1]. 
Moreover, by Proposition 14.11 we deduce that p2 = P2,o; hence p2 = 
Pi = P2,o- This solution respects all the properties of Definition 13.11 and 
the entropy condition ([6]). 

3. pifl < a < p2fi. By Proposition l3.lt we deduce that pi = pi o or pi > a 
and that p2 = P2,o or p2 < cr. 

If /(pi,o) = /(p2,o), then, by Proposition 14.11 the only possibility is 
that pi = pifi and p2 = p2,o- 

If /(pi,o) > /(P2,o)) then, by Proposition 14.11 the only possibility is 
that pi = P2 = P2,o- 

Finally, if /(pi,o) < f{p2,o), then, by Proposition 14. 11 the only possibil- 
ity is that pi = P2 = Pi,o- 

In all the cases, the solution respects all the properties of Definition l3.1l 
and the entropy condition 

4. P2,o < cr < pi^o- By Proposition 13. 11 we deduce that pi> a and p2 < cr. 
By Proposition 14.11 the only possibility is that pi = P2 = o"- The 
solution respects all the properties of Definition 13.11 and the entropy 
condition 

The proof is completed. □ 

Remark 3 In \TB^ . the authors described all the Riemann solvers, with suit- 
able properties, for nodes J with n = m = 1. The unique Riemann solver 
TZS satisfying (El) corresponds to the Riemann solver generated by the set 
X = {/(cr)} and described in Section 3.1 of J7^. 

Remark 4 One can try to generalize the entropy condition (El), at least 
for nodes with n = m = 1, to the case of fluxes depending on the arcs. 
Unfortunately this is not a trivial problem. Consider indeed the following 
example. Let /i : [0, 1] ^ M, /2 : [0, 1] — >• R be two fluxes satisfying {T) and 
assume that: 

1. fi is the flux in the arc Ii; 

2. f2 is the flux in the arc I2; 

3. cr = I the point of maximum for both fi and f2; 
4- flip) < /2(p) for every p g]0, 1[. 
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Figure 1: The situation in the example of Remark |H 

Choose < p2 < pi < ^ such that /i(pi) = /2(p2) and take k G [p2,Pi]; see 
FigureUi Then, the entropy condition ([6]) becomes 

/l(Pl)-/l(^)>/2(^)-/2(p2), 

which is equivalent to fi{k) + f2{k) < /i(pi) + /2(p2)- The last inequality 
does not hold for k = pi and for all k G [p2, pi] near to pi. 

4.2 Nodes with n = m = 2 

Consider a Riemann solver TZS for a node J with two incoming and two 
outgoing arcs. In this subsection, we assume that {pi, p2, Ps, Pa) denotes an 
equilibrium for TZS. Recall that the equilibrium must satisfy /(pi) + f{p2) = 
fips) + f{p4,)- By symmetry, we may assume also that 

(HI) pi < p2 and ps < p^. 

The results of this subsection are summarized in Table [H 

Proposition 4.2 Assume (HI) and that every pi (I G {1,2,3,4}) is a good 
datum. 

1. If TZS satisfies the entropy condition (El), then pi = p2 = Pz = Pi = cr- 

2. If pi = p2 = p3 = pi = a, then J^{pi, p2, ps, p4, k) = 0, for every 
k G [0,1]. 

Proof. Since all the data are good, then ps < p4 < cr < pi < p2. 
If /c G [p3,p4], then the entropy condition (El) becomes 

/(pi) + /(p2)-2/(A;)>/(p4)-/(p3), 
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which is equivalent to f{k) < /(ps). This imphes that f{p4,) = fips) and so 

P3 = Pi- 

If A; G [pi,p2], then in the same way we deduce that pi = p2- 

Finally, if /c G [p4, pi], then ([6]), coupled with the previous results, becomes 

2/(pi)-2/(A;)>2/(fc)-2/(p4), 

which is equivalent to f{k) < /(pi). Therefore pi = a and the conclusion 
follows. □ 



Proposition 4.3 Assume (HI) and that the equilibrium (pi, p2, ps, P4) for 

TZS is composed by three good data and one bad datum. 

1. Assume that the bad datum is in an incoming arc, say pi < a. 

If TZS satisfies (El), then P2 = <y and both p^ and p^^ belong to [pi,o"]. 
If P2 = <^ and both p3 and p^ belong to [pi, a], then T (pi, p2, Ps, P4, A;) > 
for every k G [0, 1]. 

2. Assume that the had datum is in an outgoing arc, say p4 > a. 

If TZS satisfies (El), then ps = a and both pi and p2 belong to [cr, P4]. 
If Ps = a and both pi and p2 belong to [a, P4], then T (pi, p2, Ps, P4, ^) > 
for every k G [0, 1]. 

Proof. First assume that the bad datum is in an incoming arc and the Rie- 
mann solver satisfies the entropy condition (El). Without loss of generality, 
suppose that pi < ex, p2 > a and Pa < P4 < o. We have three possibilities. 

(a) pi < Ps < P4- 

(b) Ps < Pi < P4- 

(c) Ps < P4 < Pi- 

Consider the case (a). If A; G [pi,p3], then the entropy condition (El) 
becomes 

/(p2)-/(pi)>/(ps) + /(p4)-2/(A;), 

equivalent to /(/c) > /(pi), which is true. 

If G [p4,P2], then the entropy condition (El) becomes 

/(p2)-/(pi)>2/(fc)-/(p3)-/(p4), 
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equivalent to /(P2) > f{k), which imphes that p2 = o. 
If /c G [p3,P4], then the entropy condition (El) reads 

/(P2)-/(P1)>/(P4)-/(P3), 

equivalent to /(cr) > /(P4), which is true. 

Consider the case (b). If A; G [p3,pi], then the entropy condition (El) 
reads 

/(pi) + /(p2)-2/(A;)>/(p4)-/(p3), 

equivalent to /(P3) > This implies that p3 = pi and so we are in the 

case (a). 

Consider the case (c). If A; G [p3,P4], then the entropy condition (El) 
becomes 

/(pi) + /(p2)-2/(A;)>/(p4)-/(p3), 

equivalent to /(P3) > fik)- This implies that p3 = P4. 
If A; G [p4,pi], then the entropy condition (El) reads 

/(pi) + /(p2)-2/(fc)>2/(A:)-2/(p4), 

i.e. /(P4) > /(A;). This implies that p4 = pi and so we have a contradiction 
since, by case (a), pi = pg = p4 < a = p2 and so /(pi)+/(p2) 7^ /(P3)+/(P4)- 
The second statement in the case the bad datum is in an incoming arc 
easily follows. 

Assume now that the bad datum is in an outgoing arc and that the Rie- 
mann solver satisfies the entropy condition (El). Without loss of generality, 
suppose that Ps < cr, > a and a < pi < P2- We have three possibilities. 

(a) pi < P2 < P4- 

(b) pi < P4 < P2- 

(c) P4 < Pi < P2. 

Consider the case (a). If A; G [p3,pi], then the entropy condition (El) 
becomes 

/(pi)+/(p2)-2/(A;)>/(p4)-/(p3), 

i.e. /(ps) > f{k). This implies that p3 = cr. 
If A; G [p2,P4], then ([6]) becomes 

2/(fc)-/(pi)-/(p2)>/(p4)-/(p3), 
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equivalent to f{k) > f{p4), which is true. 
If k E [pi,p2], then becomes 

/(P2)-/(P1)>/(P4)-/(P3), 

equivalent to /(a) = /(/O3) > f{pi), which is true. 

Consider the case (b). If A; G [p4, P2], then the entropy condition (El) 
reads 

/(p2)-/(pi)>2/(fc)-/(p4)-/(p3), 

which is equivalent to /(P2) > /(^)- Thus we deduce that p2 = pi and so we 
are in the case (a). 

Consider the case (c). If A; G [pi,p2], then the entropy condition (El) 
reads 

/(P2)-/(P1)>2/(A;)-/(P4)-/(P3), 

equivalent to /(P2) > f{k). This implies that pi = p2. 
If A; G [psjpJ, then ^ reads 

/(Pl)+/(p2)-2/(fc)>/(p4)-/(P3), 

i.e. /(ps) > f{k) and so pa = a. Therefore pi = P2 = P3 = P4 = c, which is 
a contradiction. 

The second statement of the item 2 of the Proposition easily follows. The 
proof is finished. □ 



Proposition 4.4 Assume (HI) and that the equilibrium {pi, P2, Ps, Pa) for 
TZS is composed by two good and two bad data. 

1. Assume that p2 < cr, i.e. the bad data are both in the incoming arcs. 
If the Riemann solver TZS satisfies the entropy condition (El), then 
Pi < P3 < p4 < p2. 

U Pi < Pi, < Pi < P2 < (y, then J^(pi,p2,p3,p4, A;) > for every 
k G [0,1]. 

2. Assume that ps > a, i.e. the bad data are in the outgoing arcs. 

If the Riemann solver TZS satisfies the entropy condition (El), then 

P3 < Pi < P2 < P4- 

If cr < Ps < Pi < p2 < Pa, then {pi, p2, ps, Pa, k) > for every 
k G [0,1]. 
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3. Assume that pi < a < p^, i.e. the bad data are in the arcs Ii and I4. 
If the Riemann solver TZS satisfies the entropy condition (El), then 
Pi < P3 < cr < p2 < Pa- 

If Pi < P3 < (T < P2 < Pa, then J^" {pi, p2, ps, Pa, k) > for every 
k e [0,1]. 

Proof. Assume that P2 < cr and that the Riemann solver satisfies the 
entropy condition (El). Since there are exactly two bad data, then p4 < cr. 
The conservation of mass at J implies that we have the following possibilities. 

(a) pl< P3<PA< P2- 

(b) P3 < Pi < P2 < Pa- 

Consider the case (a). If A; G [pi,p3], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi), which is true. 
If /c G [p3,P4], then (E]) becomes 

/(P2)-/(P1)>/(P4)-/(P3), 

which clearly holds. 

If A; G [p4,P2], then the entropy condition (El) reads 

/(p2)-/(pi)>2/(fc)-/(p3)-/(p4), 

equivalent to /(P2) > f{k), which is true. 

Consider the case (b). If A; G [p3,pi], then the entropy condition (El) 
reads 

/(Pl) + /(P2)-2/(fc)>/(p4)-/(P3), 

equivalent to /(P3) > f{k). This implies that pi = p-^ and consequently 
P2 = Pa- The second statement of the item 1 of the Proposition easily follows. 

Assume now that ps > a and the Riemann solver satisfies the entropy 
condition (El). Consequently pi > cr. Since /(pi) + /(P2) = /(ps) + /(P4), 
we have the following possibilities. 

(a) P3 < Pi < P2 < Pa- 

(b) pl< p3< PA< P2- 
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Consider the case (a). If A; G [p3,pi], then the entropy condition (El) 
reads 

/(Pi) + /(p2)-2/(A;)>/(p4)-/(p3), 

equivalent to /(ps) > f{k), which is true. 

If A; G [pi,p2], then the entropy condition (El) reads 

/(P2)-/(P1)>/(P4)-/(P3), 

which clearly holds. 

If A; G [p2,P4], then 1^ reads 

2/(fc)-/(pi)-/(p2)>/(p4)-/(P3), 

equivalent to f{k) > f{pi), which is true. 

Consider the case (b). If /c G [pi,p3], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi). This implies pi = ps and so p2 = P4. The second 
statement of the item 2 of the Proposition easily follows. 

Assume now pi < a < p4, i.e. the bad data are in the arcs Ii and 1^, and 
that the Riemann solver satisfies the entropy condition (El). We have the 
following possibilities. 

(a) Pi < P3 < 0- < P2 < P4- 

(b) P3 < pi < (7 < P4 < P2. 

Consider the case (a). If A; G [pi,p3], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi), which is true. 
If A; G [p3,p2], then (E]) becomes 

/(P2)-/(P1)>/(P4)-/(P3), 

equivalent to /(P3) > /(pi), which is true. 

If A; G [p2,P4], then the entropy condition (El) becomes 

2/(A:)-/(pi)-/(p2)>/(p4)-/(p3), 
equivalent to f{k) > f{p4), which is true. 
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Consider the case (b). If /c G [p3,pi], then the entropy condition (El) 
reads 

/(Pi)+/(p2)-2/(A;)>/(p4)-/(p3), 

equivalent to /(ps) > f{k). This imphes pi = ps and so p2 = P4. The second 
statement of the item 3 of the Proposition easily follows. 

The proof is finished. □ 



Proposition 4.5 Assume (HI) and that the equilibrium (pi, P2, Ps, P4) for 

TZS is composed by three bad data and one good datum. 

1. Assume that p2 > cr, i.e. the good datum is in an incoming arc. 

If the Riemann solver satisfies the entropy condition (El), then pi < a, 
P3 > (T, p2< P4 and /(pi) < max{/(p2),/(p3)}- 

If pi < a, Ps > a, p2 < P4 and /(pi) < max {/(P2), /(ps)}, then 
J'(pi,P2,P3,p4, k) >0 for every k G [0, 1]. 

2. Assume that p3 < a, i.e. the good datum is in an outgoing arc. 

If the Riemann solver satisfies the entropy condition (El), then p2 < o , 
Pi > cr, p3 > pi and f{p^) < max {/(pa), /(ps)}. // P2 < (x, P4 > a, 
p3 > Pi and f{pi) < max{/(p2),/(p3)}, then J=' {pi, p2, p3, Pa, k) > 
for every k G [0, 1]. 

Proof. Assume first that p2 > o" and that the Riemann solver satisfies the 
entropy condition (El). We easily deduce that pi < <J < ps < pi- We have 
the following possibilities. 

(a) pi < o- < P2 < P3 < p4. 

(b) pi < o- < P3 < P2 < Pi- 

(c) Pl < 0- < P3 < P4 < P2. 

Consider the case (a). If A; G [pi,p2], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi). This implies that /(P2) > fi.Pi)- 
If /c G [P25P3], then ^ becomes 

2/(fc) - f{pi) - fip^) > /(P3) + /(P4) - 2f{k), 
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equivalent to 2f{k) > /(ps) + fiPi), which is true. 
If k E [p3,P4], then the entropy condition (El) reads 

2/(A:)-/(pi)-/(p2)>/(p4)-/(P3), 

equivalent to f{k) > f{pi), which is true. 

Consider the case (b). li k E [pi,p3], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi). This implies that /(P3) > /(pi)- 
If A; G [p3,P2], then ([6]) becomes 

/(P2)-/(P1)>/(P4)-/(P3), 

equivalent to f{p3)>f{pi)- 

If A; G [p2,P4], then the entropy condition (El) reads 

2/(A;)-/(pi)-/(p2)>/(p4)-/(P3), 

equivalent to f{k) > /(P4), which is true. 

Consider the case (c). If A; G [p4,P2], then the entropy condition (El) 
reads 

/(P2)-/(P1)>2/(A;)-/(P3)-/(P4), 

equivalent to /(P2) > f{k). This implies that p2 = P4 and so we are in 
the case (b). The second statement in the item 1 of the Proposition easily 
follows. 

Assume now that ps < a and that the Riemann solver satisfies the entropy 
condition (El). We easily deduce that pi < P2 < cr < P4. We have the 
following possibilities. 

(a) pi < P2 < P3 < o" < P4. 

(b) Pi < p3 < p2 < o- < p4. 

(c) P3 < pi < P2 < cr < P4. 

Consider the case (a). If A; G [pi,P2]! then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi), which is true. 
If A; G [p2,P3], then ([6]) becomes 

2f{k) - /(pi) - /(p2) > /(P3) + /(P4) - 2/(A;), 
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equivalent to 2f{k) > f{pi) + f{p2), which is true. 
If k E [p3,P4], then the entropy condition (El) reads 

2/(A;)-/(pi)-/(p2)>/(P4)-/(p3), 

equivalent to f{k) > /(P4). This implies that /(ps) > f{pA)- 

Consider the case (b). If A; G [pi,p3], then the entropy condition (El) 
reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to f{k) > /(pi), which is true. 
If /c e [p3,P2], then dH]) becomes 

/(P2)-/(P1)>/(P4)-/(P3), 

equivalent to /(P3) > /(pi), which is true. 

If A; G [p2,P4], then the entropy condition (El) reads 

2/(fc)-/(pi)-/(p2)>/(p4)-/(p3), 

equivalent to f{k) > fip^j- This implies that /(P2) > fiPi)- 

Consider the case (c). If A; G [p3,pi], then the entropy condition (El) 
reads 

/(Pl)+/(P2)-2/(fc)>/(p4)-/(p3), 

equivalent to /(P3) > /(A;). This implies that pi = p3 and so we are in 
the case (b). The second statement in the item 2 of the Proposition easily 
follows. 

The proof is finished. □ 

Proposition 4.6 Assume (HI) and that the equilibrium [pi, p2, Ps, P4) for 

TZS is composed by four bad data. If the Riemann solver satisfies the entropy 
condition (El), then pi < P2 < cr < Ps < P4- Moreover, if pi < P2 < cr < 
P3 < P4, then T (pi, P2, P3, p4, A;) > for every k G [0, 1]. 

Proof. It is sufficient to check the entropy condition (El). If A; G [pi,P2], 
then the entropy condition (El) reads 

/(p2)-/(pi)>/(p3) + /(p4)-2/(A;), 

equivalent to /(A;) > /(pi), which is true. 
If A; G [P25P3], then (E]) becomes 

2f{k) - /(pi) - /(p2) > /(P3) + /(P4) - 2/(A;), 
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Bad data 


admissible configurations 





Pi = P2 = P3 = P4 = cr 


1 


Pi < P3 < P4 < O- = P2, Pi < cr 


P3 = 0- < Pi < P2 < P4, P4 > cr 


2 


Pl < P3 < P4 < P2 < 0- 


0- < P3 < pi < P2 < P4 


Pi < P3 < 0- < P2 < P4, Pl < 0- < P4 


3 


pi < a < p3 < p4, cr < p2 < P4, /(pi) < max {/(P2), /(ps)} 


Pl < P2 < (T < P4, Pl < P3 < (T, /(p4) < max{/(p2), /(ps)} 


4 


Pl < P2 < 0- < P3 < P4 



Table 1: All the possible configurations for an equilibrium (pi, P25 Ps, P4) of 
a TZS satisfying the entropy condition (El). By symmetry, we assume that 
Pl < P2 and p3 < p4, i.e. (HI) holds. 



equivalent to 2/(A;) > /(ps) + /(P4), which is true. If A; G [p3,P4], then the 
entropy condition (El) reads 

2/(A;)-/(pi)-/(p2)>/(P4)-/(p3), 
equivalent to f{k) > f{p4), which is true. This concludes the proof. □ 

Remark 5 Note that there exist Riemann solvers satisfying the consistency 
condition and the entropy condition (El). Here we construct a Riemann 
solver TZS with such properties. 

Consider an initial condition (pi,o, P2,05 Ps.o, P4,o)- Denote with {pi, P2, Ps, P4) 
the image of the initial condition through TZS, i.e. 

(Pl,P2,P3,P4) = 7^>S(pi,o,P2,0,P3,0,p4,o) 

// h is the number of bad initial data, then we define TZS according to the 
following possibilities. 

h = 0. We put Pl = P2 = Ps = p4 = a. By Proposition ^.^ this provides an 
entropy admissible equilibrium. Moreover 

TZS (7^5(pi_o,P2,0,Ps,0,P4,o)) = ^'5(pi,o,P2,0,P3,0,P4,o)- 

h = 1. Let I G {1,2,3,4} be such that pjQ is a bad datum. We have two 
possibilities: I <2 orl>3. 
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Assume first J < 2. We put pj = pj^ and pi = a for I G {1, 2}, / ^ I. 
Moreover we define ps = pi and P4 = P2- 

Assume now I > 3. We put pj = pjQ and pi = a for I G {3,4}, / 7^ /. 
Moreover we define pi = pa and p2 = P4- 

By Proposition \4.3\ these solutions provide entropy admissible equilib- 
ria. Moreover 

nS (JZSipifi, P2fi, P3fi, Pifi)) = ^'5(pi,o,P2,0,P3,0,P4,o)- 

h = 2. Let li, I2 E {I, 2, 3, 4}, Zi 7^ I2, be such that pi^^ and pi^^ are bad data. 
We have three different possibilities. 

Assume first that li,l2 G {1,2}. In this case we put pi^ = pi^^, pi^ = 
Pi2,o, P3 = Pi and Pi = P2- 

Assume now that h^h G {3,4}. In this case we put pi^ = pi^^, pi^ = 
Phfi, Pi = h (^nd p2 = Pa- 

Consider finally the last case. For simplicity suppose that /i = 1 and 
/2 = 4. We define pi^ = pi^^, pi^ = pi^^, = Ph (^^d pg = p;,. By 
Proposition \4-4\ these solutions provide entropy admissible equilibria. 



Moreover 

lis (7^5(pi,o,P2,0,P3,0,P4,o)) = ^'5(pi,o,P2,0,P3,0,P4,o)- 

h = 3. Let I G {1,2,3,4} be such that p^-g is a good datum. We have two 
possibilities: I <2 or I > 3. 
Assume first I < 2; say I = 2 for simplicity. 

If fipsfi) + fipAfl) - fipifl) e [min{/(p3,o_),/(p4,o)},/(o^)], then we 
put pi = pi Q for every I G {1,2,3,4}, / 7^ / and pf G [a, 1] such that 

f{p2) = /(P3,o) + /(P4,o) - /(Pi,o)- 

If fipsfi) + /(P4,o) - /(Pi,o) > /(o-) and /(p3,o) > /(P4,o); then pi = 
P3 = Pi,o and p2 = Pa = P4,o- 

If f{p3,o) + /(P4,o) - /(pi,o) > /(o-) and /(ps.o) < /(p4,o), then p2 = 
P3 = P3,o and pi= Pa = pi,o- 

If fip3,o) + fip4fl)-fipi,o) < min{/(p3,o),/(P4,o)}, then p2 = Pa = P3,o 
and pi = Ps = P4,o- 

Assume now I > 3; say I = 3 for simplicity. 

If fipifi) + /(P2,o) - /(P4,o) e [min{/(pi,o),/(p2,o)},/(c^)], then we 
put pi = pifi for every I G {1,2,3,4}, / 7^ / and pj G [0, cr] such that 

fips) = /(Pi,o) + /(P2,o) - /(P4,o)- 

If fipifl) + /(P2,o) - /(P4,o) > /(o-) and /(pi,o) > /(P2,o); then pi = 
P4 = p4,o and P2 = P3 = P2,o ■ 

If fipifi) + /(P2,o) - /(P4,o) > /(o") and /(p2,o) > /(pi,o), then p2 = 
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Pa = p4,o and pi = = Pi,o- 

U f{pi,o) + f{p2,o)-f{pA,o) < min{/(pi,o),/(p2,o)}, then p2 = Pa = P2,o 
and pi = P3 = Pi,o- 

By Propositions \4-4\ and \4-5\ these solutions provide entropy admissible 



equilibria. Moreover 

TZS (7^5(pi,o,P2,0,P3,0,P4,o)) = ^'5(pi,o,P2,0,P3,0,P4,o)- 

h = A. We have some different cases. Assume first that /(pi,o) + /(P2,o) = 
/(P3,o) + /(p4,o)- We put pi = pi,o, P2 = p2,o, P3 = p3fi and pi = p4,o- 
Assume now that /(pi,o) + /(p2,o) < /(P3,o) + /(p4,o)- For simplicity 
suppose that /(pi,o) < /(P2,o) and /(p3,o) > /(P4,o)- 
V f{pi,o) > f{p2fi), then we put p^ = p^ = p^ o and p2 = p^ = p2,o- 
U /(P4,o) < /(P2,o); then we put pi = pi^o, pa = p2,o, Pa = P4,o and 
p3 e [0,(t] snc/i i/iai /(pg) = /(pi) + /(p2) - fipi). 

Assume finally that /(pi,o) + /(P2,o) > /(P3,o) + /(P4,o)- For simplicity 

suppose that /(pi,o) < /(P2,o) and /(p3,o) > /(P4,o)- 

If /(pi,o) > /(P3,o); then we put p^ = p^ = p^ ^ and p2 = Pa = Pa,o- 

If /(pi,o) < /(P3,o), then we put pi = pi,o, P3 = P3,o, P4 = P4,o and 

P2 e [cr, 1] such that /(pa) = /(P3) + /(P4) - /(pi)- 

By Propositions \4-4\ \4-5\ and \4.6\ these solutions provide entropy ad- 
missible equilibria. Moreover 

nS (7^5(pi,o,P2,0,P3,0,P4,o)) = ^'5(pi,o,P2,0,P3,0,P4,o)- 

Remark 6 Another example of Riemann solver satisfying the entropy condi- 
tion (El) for a node with two incoming and two outgoing arcs is a particular 
case of the Riemann solver IZS2, defined in Section \5T^ see Proposition \5.5[ 
The Riemann solver TZS, constructed in Remark\^ differs from the Rie- 
mann solver TZS2 ■ The key difference is that a permutation of initial data 
in incoming (resp. outgoing) arcs influences the solution in outgoing (resp. 
incoming) arcs in the case of TZS, but not in the case ofTZS2. 
Consider the following example. Let f{p) = 4p(l — p) be the flux. Assume 
that (|,|,;|,|) and (|, |, |, |) are two initial conditions. In both cases, we 
have only one bad datum and so, using the notation of Remark{5l h = 1. 
Hence we deduce 



TZS 

and 



1 3 1 1\ /I 1 1 1 



4'4'4'4/ V4'2'4'2 



4' 4' 4' 4 / V 2' 4' 2' 4 ' ' 
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while 

1 3 1 1\ _^ /3 1 1 1 



' 4'4'4'4 j V4'4'4'4 " 



see Section 



5 Examples 

This Section deals with some examples of Riemann solvers, introduced in 
literature for describing car and data traffic. For each of them, we analyze 
the entropy conditions (El) and (E2). First we need some notation. 
Consider the set 

{0 < aji < 1 Wi,j, ^ 
A = {a,,} : "f™^ _ivz \ (^4) 

j=n+l,...,n+m — 1 Vt 

j=n+l ) 

Let {ei, . . . , e„} be the canonical basis of M". For every i = 1, . . . , ra, we 
denote Hi = {ci}^. If A & A, then we write, for every j = n + 1, . . . ,n + m, 
ttj = {ttji, . . . ,ajn) € M" and Hj = {aj}^. Let /C be the set of indices 
k = [ki, ki), 1 < i < n — 1, such that 0<ki<k2<---<ki<n + m 
and for every k G /C define 

e 

h=l 

Writing 1 = (1,...,1)g]R" and following [H] we define the set 

'n:= {Ae A:l ^ for every k G /C} . (15) 

Notice that, if n > m, then 91 = 0. The matrices of will give rise of a 
unique solution to the Riemann problem at J. 
For later use, define the set 

5.1 Riemann Solver TZSi 

In this subsection, we consider the Riemann solver introduced for car traffic 
in [6]. The construction can be done in the following way. 
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1. Fix a matrix A and consider the closed, convex and not empty set 



n+m 



I i=l j=n+l ) 

2. Find the point (71, . . . , 7^) G Vt which maximizes the function 

E(7i,...,70 =7i + --- + 7„, (18) 

and define (7n+i, . . . , 7n+m)^ := A - (71, ... , 7„)^. Since A G DT, then 
(71, . . .,%) is unique. 

3. For every i G {1, . . . , n}, define pi either by pi^ if f{pifi) = 7j, or by the 
solution to /(p) = 7j such that pi > a. For every j G {^+1, • • • , n+m}, 
define pj either by pj^ if f{pjfl) = 7j, or by the solution to /(p) = 7^ 
such that Pj < a. Finally, defi^ne TZSi : [0, 1]"+"^ [0, 1]'^+'" by 

,0; • • • ) Pn+m,0 ). (19) 

The following result holds. 

Lemma 5.1 The function defined in / figj) satisfies the consistency condition, 
in the sense of Definition \3. 5l 

For a proof, see P [IZ] . We show that this Riemann solver does not satisfy 
neither the entropy condition (El) nor (E2). 

Proposition 5.1 The Riemann solver TZSi does not satisfy the entropy con- 
dition (E2) in the sense of Definition \3. 6\ and, consequently, does not satisfy 
the entropy condition (El) in the sense of Definition \3. 5i 

Proof. Consider a node with 2 incoming and 2 outgoing arcs, the fiux 
function /(p) = 4p(l — p), a matrix 



A 



1 1 
3 2 

2 1 

3 2 



and the initial conditions pi^ = |, p2,o = |, P3,o = Pi,o — Jq- 

this case the set Q in ( |T7|) is 



(71,72) e [0,1] X 



0.^ 

' 16 



7i 72 15 271 72 

:0< — + — <— , 0< — + — <1 
-3 2 - 32 3 2 - 
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72 



ft 



7i 



Figure 2: The set fl of Proposition 15. 1[ 



see Figure 121 Therefore we deduce that 7i = 72 = If, 73 = If, 74 = 
pi = a, p2 > (J, P3 = psfl and p^ < a. The entropy condition ([7]) in this case 
becomes 

/(P2) - /(a) > /(p3) - f{o) + /(a) - /(p4), 
which is equivalent to 

IS 1 ^ 77 1 Q 

0</fe)-/M-/fe) + /fe) = j3-l-3^ + gg = -5. 



This concludes the proof. 



□ 



The maximization of the function E over which defines the Riemann 
solver 7^i5i, is, however, in connection with the maximization of the entropy 
T . In order to explain this fact, let us introduce some notations. 
Given f2 in (ITTIl . define 



(/(P 



1^ 



,/(Pn)) G 



n+m 



{pi, . . . , pn+m) ^ Yl^l 



1=1 



( f{Pn+l) \ 




( /(Pi) 


\ 




= A- 






\ f{pn+m) J 




V /(p«) 


/ 



> (20) 



and the functional 



g : $ 

(Pl5 • • • 5 Pn+r. 



•^(Pl, • • • ,Pn+m,0-) 



(21) 



which is the restriction of on $ x {a}. Note that the set $ consists in all 
the possible solutions at J satisfying Definition 13. II and the distribution rule, 
determined by the matrix A G Dl. It is easy to see that there exists a one to 
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one correspondence between and 

For every H C {1, . . . , n + m} of cardinality h, with < /i < n — 1, define 



(7n+l, • • • , Jn+mV = A-{ji, -fnf 



nn = I (71, . . . , 7n) e n ■ ■ ■ o^^i ^ Yrl^' 



i=l 



■ji = max fli if Z e 
ji < max if I ^ H, 



(22) 



and 



(/(pi),...,/(p„)) G (^w, 



(Pl,---,P: 



n+m 

■n+m) e n • 
Z=l 



/ fiPn+l) \ 




( /(Pi) \ 




= A- 




\ f{Pn+m) J 




^ fipn) J 



} (23) 



Notice that D,n and depend on the initial condition (pi,o, ■ ■ ■ , Pn+m,o) and 
on the matrix A G There is a one to one correspondence between Q-^ and 
given by the one-to-one function 

— > 

{Pl, ■ ■ ■ , Pn+m) I ' {f(Pl),---,f(Pn))- 

Moreover, if fl-^ ^ 0, then has, at most, topological dimension n — h. 
The following proposition holds. 

Proposition 5.2 Let Ti C {1, . . . ,n + m} &e a sei o/ cardinality h, with 
< h < n — 1 and suppose that fl-^ ^ 0. The functional Q, restricted to 
is given by 

G{pU...:Pn+m)^ E [fiPl)-fi'^)]+J2^fi'^)-fiPl)]- ^ 

ie{i,...,n+m}\n leH 

Proof. Fix (pi . . . , Pn+m) e and I e {1, . . . , n -I- m}. We have some 
different possibilities. 

1. I < n and I E H. In this case the term sgn(p; — a) {f{pi) — /(cr)) 
becomes /(cr) - f(pi). 

2. I < n and I ^ 7i. In this case the term sgn(p; — a) (/(p;) — /(o")) 
becomes f{pi)- f{a). 

3. / > n + 1 and I eH. In this case the term — sgn(p; — a) {f{pi) — /(c)) 
becomes /(a) - f{pi). 
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4. I > n + 1 and I ^7i. In this case the term — sgn(pi — a) {f{pi) — /(o")) 
becomes fipi)- /(a). 

Therefore the proof is finished. □ 

Corollary 5.1 Let Ti C {1, . . . , n + m} he a set of cardinality h, with < 
h < n — 1 and suppose that Q-^ ^ 0. The problem of maximizing Q on the 
set is equivalent to the problem of maximizing the function E, defined 
in on the set fin- 

Proof. Notice that, by Proposition 15.21 the function Q on the set 
coincides with 

l&{l,...,n+m}\H 

where C is a constant, depending on Ti. and on the initial conditions. Indeed, 
if / G then pi is completely determined by the initial condition pi^o- More 
precisely, pi is equal to pi^ when pi^o is a bad datum, while pi is equal to a 
in the other case. Therefore, if (pi, . . . , Pn+m) ^ then we deduce that 

g{pi, pn+m) = Yl f^P^^ + 5Z f^Pi^ + 

j6{l,...,n}\'K je{n+l,...,n+m}\n 

= E f^p^^+ E fip^)+Ci 

ie{l,...,n}\'H je{n+l,...,n+m} 

= E f(p^^+ E 

ie{l,...,n}\W iG{l,...,n} 

= 2 E f(P^) + C2, 

ie{l,...,n}\W 

where Ci and C2 are constants. Finally note that the function E, restricted 
on Q-j-i, is given by 

E(7i,...,7„)= E 7. + C2-C1. 

je{l,...,n}\W 

This completes the proof. □ 

Remark 7 Note that the set $ is, in general, disconnected, while the set fl 
is convex and so connected. The function Q, defined in [21\) . i.e. the entropy 
function restricted on $ x {a}, is continuous, since it has not jumps in each 
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connected component of^. Since there is a bijection between the sets and 
$, then we can consider the entropy function on fl. More precisely, define 
the function 

T -.n — > $ 

(7l,---,7n) ' ' (Pl, • • • ,Pn+m), 

satisfying f{pi) = 7i for every i G {1, . . . , n}, and consider the map Q oT : 
f2 — i> R. This map, in general, is discontinuous, since it can have jumps 
at every point (71, • • • ,7n) £ ^ ^^^^ ^1 '^2 different subsets of 
{1, . . . , n + m} of cardinalities less than or equal to n — 1. 



5.2 Riemann Solver IZS2 

In this subsection, we consider the Riemann solver, introduced in [13] for 
data networks; see also [17] . The construction can be done in the following 
way. 

1. Fix E Q and define 

n n+m 

Tine = sup ^i, Tout= ^ sup^-, 

i=l j=n+l 

then the maximal possible through-fiow at the crossing is 

r = min {Tine, Tout} ■ 

2. Introduce the closed, convex and not empty sets 

{n n 
(7i,---,7n)en^- 5Z^^ = r 
i=l i=l 

{n+m n+m ^ 

(7n+l, • • • , 7n+m) € JJ : ^ = \ ■ 

j=n+l j=n+l J 

3. Denote with (71, . . . ,7^) the orthogonal projection on the convex set 
/ of the point {T9i, . . . , T9n) and with (7„+i, . . . , 7n+m) the orthogonal 
projection on the convex set J of the point {T9n+i, ■ ■ ■ , T9n+m)- 

4. For every i G {1, . . . , n}, define pi either by pi^ if f{pifl) = 7j, or by the 
solution to /(p) = 7j such that > a. For every j G {n+1, . . . , n+m}, 
define pj either by pj^ if f{pjfl) = 7j, or by the solution to /(p) = 7^ 
such that Pj < a. Finally, defi^ne : [0, 1]"+™ ^ [0, 1]'^+"^ by 

) = (Pl,... ) Pn, Pn+1, ■ ■ ■ , Pn+m ) • (25) 
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The following result holds. 
Lemma 5.2 The function defined in f25\} satisfies the consistency condition 

nS 2(113 2{pi,0, Pn+mfl)) = ^'52(pi,0, • • • , Pn+mfl) (26) 

for every (pi,o, . . . , pn+m,o) e [0, 1]"+"". 

For a proof, see [H] . We prove now that the Riemann solver TZS2 satisfies 
the entropy condition (E2). 

Proposition 5.3 Assume n = m and consider a node J with n incoming 
roads and m outgoing roads. The Riemann solver TZS2 satisfies the entropy 
condition (E2) in the sense of Definition \3. 6i 

Proof. Fix an initial condition (pi,o, • • • , Pn+m,o) and define (pi, . . . , p„+m) = 
7^.52 (pi,o 5 • • • 5 Pn+m,o)- Wc have two different cases. 

Fine < Tout- In this situation, we deduce that pi < a for every z G {1, . . . , n}. 
Thus the entropy reads 

n n+m 

T{pi, pn+m, (t) = n/((T)-^ f{pi)-Y^ sgn(pj— a) (/(p^) - /(a)) . 

i=l j=n+l 

For every j E {n+1, . . . ,n + m}, the term — sgn(pj — cr) {f{pj) — f{cr)) 
can be minorized by f{pj) — f{cr) and so 

n n+m 

J^(pi,...,p„+^,a) > - ^/(pi) + ^ (/(pj) - /(a)) 

1=1 j=n+l 

= (n — m)f{a) = 0. 

Tine > Tout- In this situation, we deduce that pj > a for every j G {n + 
1, . . . ,n + m}. Thus the entropy reads 

n n+m 

T{pi,...,pn+m,(y) =^sgn{pi-a){f{pi)-f{a))+mf{a)-Y^ f{pj). 

i=l j=n+l 

For every i G {1, . . . , n}, the term sgn(pj — a) {f{pi) — /(o")) can be 
minorized by /(pj) — /(a) and so 



n+m 



J^{pi,...,Pn+m,(^) > ^ifipi) - fia)) + mf{a) - ^ /(p 

i=l j=n+l 

= (m — n)f{a) = 0. 
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The proof is finished. □ 

In general, the Riemann solver TZS2 does not satisfy the entropy condition 
(El) even in the case n = m, as the next Proposition shows. 

Proposition 5.4 The Riemann solver TZS2 does not satisfy the entropy con- 
dition (El) in the sense of Definition \3.5[ 

Proof. Consider a node with 2 incoming and 2 outgoing arcs, the flux 
function /(p) = 4p(l — p), = (i, i, ^, ^) and the equilibrium configuration 

(^i' hl~ 1 ~ 471) ■ ^^^^ equation becomes 
-sgn (l--l--k] (l-fik)) >0 



2 aV2 

for every G [0, 1]. If = j, then the previous inequality becomes 

8 4; VS 4; - ' 
which is clearly false. □ 

Indeed, in some special situation, namely for nodes with 2 incoming and 
2 outgoing arcs and = (|, |, ^, |), the Riemann solver TZS2 satisfies the 
entropy condition (El). 

Proposition 5.5 Fix a node J with two incoming and two outgoing arcs. If 
= (i, i, i, i), then the Riemann solver 7IS2 satisfies the entropy condition 
(El), in the sense of Definition \3.5\ 

Proof. Consider an equilibrium {pi, P2, Ps, Pa) for the Riemann solver 71S2 
and denote with g the number of good data. We have the following possibil- 
ities. 

(7 = 4. In this case we deduce that {pi, P2, Ps, Pi) = (i^ |? |? i) so the 
entropy condition (El) is satisfied. 
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(7 = 3. Consider only the case T = Tine, since the other case T = Tout is 
completely symmetric. Thus the bad datum is in an incoming arc and 
so we may assume that pi < cr, p2 > cr and ps < P4 < cr. Since 
6 = (|, |, |, I), then p2 = cr and Pz = pi < cr. Moreover, the fact that 
fi.pl) + /(P2) = /(ps) + /(P4) implies that 

pi < p3 = p4 < p2 = (y. 

By item 1 of Proposition 14.31 the entropy condition (El) holds. 

g = 2. Consider only the case T = Tine, since the other case T = Tout is 
completely symmetric. We have two possibilities: either the bad data 
are in the incoming arcs or one bad datum is in an incoming arc and 
the other bad datum is in an outgoing arc. 

Assume first that the bad data are in the incoming arcs. Without loss 
of generality we may assume that pi < P2 < cr and pa < p4 < a. Since 
= (i, |, |, I), then p3 = p4 and so, the fact that /(pi) + /(P2) = 
/(Ps) + fiPi) implies that 

pi< p3 = p4< p2 < cr. 

By item 1 of Proposition 14.41 the entropy condition (El) is satisfied. 

Assume now that one bad datum is in an incoming arc and the other 
bad datum is in an outgoing arc. Without loss of generality we may 
assume that pi < a < pi and ps < c < p2- Since F = Tine, then we de- 
duce that p2 = cr. Moreover 6 = (1, i, 1, implies that /(pa) > fipi) 
and so f{pi) < f{p4,), since /(pi) + /(ps) = /(ps) + fipi)- Therefore 

Pi < Ps < P2 = cr < P4 and pi < p2. 

By item 3 of Proposition 14.41 the entropy condition (El) is satisfied. 

g = I- Consider only the case F = Tine, since the other case F = Tout is 
completely symmetric. We have two possibilities: the good datum 
is in an incoming arc or in an outgoing arc. Assume first that the 
good datum is in an incoming arc. Without loss of generality, we 
may consider that pi < cr < p2 and cr < p^ < p^. Since F = Tine, 
then p2 = cr. Moreover /(pi) + /(P2) = /(Ps) + fiPi) implies that 
/(P4) > /(pi)- By item 1 of Proposition 14.51 the entropy condition 
(El) is satisfied. 

Assume now that the good datum is in an outgoing arc. Without 
loss of generality, suppose that pi < P2 < cr, p3 < cr < p4. Since 



31 



^ = (i h h ^^^^ /(Ps) > /(P4) and so /(P4) < /(ps) and pa > pi 
since f{pi) + /(P2) = /(Ps) + fiPi)- By item 2 of Proposition I4.5|, the 
entropy condition (El) is satisfied. 

= 0. In this case we have that F = Tine = ^out- Without loss of generahty, 
suppose that pi < P2 < cr < Ps < P4 and we conclude by Proposi- 
tion iJl 

The proof is finished. □ 



5.3 Riemann Solver TIS^ 

In this subsection, we consider the Riemann solver, introduced in [26] for 
crossing nodes. Consider a node J with n incoming and m = n outgoing 
arcs and fix a positive coefficient Pj, which is the maximum capacity of the 
node. The construction can be done in the following way. 

1. Fix G 0. For every i G {!,..., n}, define 

Pj = min {sup Qi, sup fii+„} . 
Then the maximal possible through-flow at J is 

n 
1=1 



2. Introduce the closed, convex and not empty set 

f n n 

I = { (71, • • • , 7n) e n[o, r.] ■Y.^^ = "^i^ Pj} 



i=l i=l 



3. Denote with (71, . . . , 7^) the orthogonal projection on the convex set / 
of the point (min{P, Tj}9i, . . . , min{P, Pj}6'„) and set (7n+i, • • • , 72n) = 

(7l,---,7n)- 

4. For every z G {1, . . . , n}, define pi either by p^ if f{pi,o) = li, or by the 
solution to /(p) = 7j such that pi > a. For every j G . . . , ra+m}, 
define pj either by pj^ if f{pj,o) = 7j, or by the solution to /(p) = ■jj 
such that Pj < a. Finally, define TZSa : [0, 1]"+™ [0, 1]"+"^ by 

) = (Pl,... ) • (27) 



32 



The following result holds. 
Lemma 5.3 The function defined in satisfies the consistency condition 

nS^inS^ipifi, Pn+m,o)) = Tl<Ss{pi^o, Pn+mfl) (28) 

for every (pi,o, . . . , pn+m,o) e [0, 1]"+"". 
For a proof, see Proposition 2.4 of 



Example 1 Consider a node J with 2 incoming arcs and 2 outgoing ones, 
6 = (|, ^, |, ^) and Tj = ||. Moreover, assume that f{p) = 4p(l — p). 
We easily see that 



{pi, P2, P3, Pa) 



11 1 /59 4 1 1/59 
5' 2 ToV T' 5' 2 ~ 10 



is an equilibrium for TZS 3 . Thus we have 

J^{pi,p2,P3,pA,(^) = (/W - /(Pi)) + (/(P2) - /(fx)) 

- (/(Ps) - /M) - (/(^) - /(P4)) 
64 

Example 2 Consider a node J with 2 incoming arcs and 2 outgoing ones, 
= (|, |, |, I) and Tj = |. Moreover, assume that f{p) = 4p(l — p). 
We easily see that 

(p-i,p-2,p-3,p-4) ={\+i\ii^i+ ^ vi' ^ + ^/i' ^ " yi 

is an equilibrium for TZS 3 . Thus we have 

■^(Pl,P2,P3,P4,f^) = (/(pi) - /((x)) + (/(pa) - /((x)) 

- (/(Ps) - /M) - (/(a) - /(p4)) 

= 2(/(p-2)-/(a)) = -| 

The following result follows by the previous examples. 

Proposition 5.6 The Riemann solver TZSs does not satisfy neither the en- 
tropy condition (El) nor the entropy condition (E2). 
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